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Fundamental quantities of QCD, such as the strong coupling and Aqqd, are studied in the frame- 
work of lattice QCD with N/ — 2 twisted mass fermions. In particular, the contact between lattice 
and continuous calculations is made by comparing the renormalized ghost-gluon vertex in MOM 
scheme with 4-loop perturbative results. A power correction is needed in order to have agree- 
ment between the two descriptions. This suggests the presence of a dimension-two (A 2 } gluon 
condensate whose value is found to be higher than in the quanched case. 
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1. Introduction 

The lattice provides a very elegant way of calculating renormalized observables. In this frame- 
work several methods are known to extract the running of the QCD coupling constant, which allows 
for the determination of the QCD scale, Aqcd and for the study of infrared properties. In the case 
of quenched world the mismatch between the perturbative running and the lattice one has revealed 
the presence of a non-null gluon condensate of dimension two that, being non-gauge invariant, has 
motivated the research of its possible implications for the gauge-invariant world. 

In this note we apply the already established methods for Nf = 2 dynamical quarks, includ- 
ing light up and down quarks. Nf = 2 + 1 + 1 lattice simulations are already being performed, 
thus a realistic lattice estimate of Ajjg directly comparable with experimental results will become 
inmediatly accesible. 

In particular here we focus on the study of the ghost-gluon vertex in the configuration of 
vanishing incoming ghost-momentum. Only in this case the ghost-gluon vertex can be related 
directly to the bare and ghost propagators, making calculations simpler. 



2. Taylor scheme 
2.1 Definitions 

In [1] was shown that the so-called Taylor scheme is the only one where the coupling can be 
cumputed from two-point Green functions, due to Taylor's theorem. We write Landau gauge gluon 
and ghost propagators as: 

(F«>)V,A) = -iu^>: (2.1) 

with A the regularisation cutoff. The renormalized dressing functions, G R and F R are defined 
through : 

G R {p\v 2 ) = limZ 3 - 1 ( J u 2 ,A)G(/,A) 

A— >o° 

F R (p 2 ,n 2 ) = limZ 3 -V,A)/V,A), (2.2) 

A— >°° 

with MOM renormalization condition 

G R {pi 2 ,^)=F R (ii\pi 2 ) = l . (2.3) 

Due to Taylor's non-renormalization theorem, the renormalized coupling defined from the ghost- 
gluon vertex with a zero incoming ghost momentum can be computed from ghost and gluon 
propagators using: 

«r(M 2 ) = = Km ^^G( A i 2 ,A 2 )F 2 ( Ai 2 ,A 2 ) ; (2.4) 

what has been called Taylor 1 scheme [1] 



'From now on, the quantities expressed in this scheme will carry the T index. 
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2.2 Perturbation theory and OPE 

The perturbative running of aj is known up to four loops [2], 




M^^fl-l^ + I^UlogW-^) + ^~-A\ (2-5) 



(M 4 \ 2j3o 2 VA 



3 



1 ( ft , 1 fft X / „, , M , , 



£r + o 7T -21og 3 (0+51og 2 (0+ log(0-l (2.6) 




with t = In ^j- and the perturbative coefficients: 

ft =j3 1 = 102- T AT / 

i 8 2 = ]8 2 -A 1 c 1 +A (c 2 -c?) 

= 3040.48 - 625.387 N f + 19.3833 N 2 - 

ft = £ 3 " 2fe + i8 i^i + A (2 c 3 - 6 c 2Cl + 4 cf ) 

= 100541 - 24423.3 + 1625.4 Nj - 27.493 #j , (2.7) 

The parameters Aqcd in two schemes can be perturbatively related at high energy. In particu- 
lar, from the T-scheme to MS this relationship reads: 



Cl 507-40N/ 
^ms = ~2f = ~192-4W f . (28) 
At 

Following the Operatore Product Expansion (OPE) program both ghost and gluon propagators 
show the appearance of a non-perturbative power correction driven by the non-gauge invariant 
dimension-two gluon condensate (see [1], [3] and referencies therein). Including power corrections 
at tree-level in ghost and gluon dressing functions, one can rewrite (2.4) as 



a r (M 2 ) = «rV)|l + 



9 gUll)(A 2 )R,4 



\x 2 A(N C - 1 



(2.9) 



where q\ ^> Aqcd is some perturbative scale and the running of the perturbative part is described 



by equation (2.5). This formula will be used for the data analysis in the next section that does 
depend on two parameters, Aqcd and (A 2 ), that will be fitted. 



3. Lattice setup and role of H(4) orbits 

The results presented here are based on the gauge field configurations generated by the Euro- 
pean Twisted Mass Collaboration (ETMC) with the tree-level improved Symanzik gauge action [4] 
and the twisted mass fermionic action [5] at maximal twist, discussed in detail in refs. [6]- [9]. 
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We preliminarly exploited 100 ETMC gauge configurations obtained for jS = 3.9 (jU = 0.0085), 
60 for j8 = 4.05 (n = 0.006) and 100 for j8 = 4.2 (/i = 0.002) simulated on 24 3 x 48 lattices, corre- 
sponding to Nf = 2 in order to compute the gauge-fixed 2-point gluon and ghost Green functions. 

For fixing Landau gauge in the lattice we minimise the functional 

Fu[g] = Re ^_ i {l-^g(x)U^x)g{x + ^ (3.1) 

respect to the gauge transform g. Ghost propagator is computed in Landau gauge as the inverse of 
the Faddeev-Popov operator, that is written as the lattice divergence, 

M{U) = ~V-D(U) (3.2) 

where the operator D acting on an arbitrary element of the Lie algebra, 77 reads: 

D(U)ri(x) = ^(U^x)ri(x + v)-ri(x)U^x) + ri(x + v)ul-ul(x)ri(x)) . (3.3) 

More details on the lattice procedure for the inversion of Faddeev-Popov operator can be found 
on [10]. 

As we intend to fit the running of a s , our interest is to have, on one hand the highest momenta 
accesible and, on the other the highest number of data points to perform the fit. When working 
at a given lattice spacing, the momentum window has to be limited due to the presence of high 
discretization errors. These lattice artifacts are due to the breaking of the rotational symmetry of 
the euclidean space-time when using an hypercubic lattice, where this symmetry is restricted to the 
discrete H(4) isometry group. These artifacts can be illustrated as the difference between the lattice 
momenta, 

p t i = ^siaap, l (3.4) 

and the continuum ones, 

Pn = J^ n = 0,l,--,N. (3.5) 

Clearly these two momenta will differ except in the limit n/N — > 0. Following what was recently 
discussed in [11] and [12], let us consider an adimensional lattice correlation function Q that de- 
pends on the lattice momentum ap^ and some mass scale aA: Q = Q(a 2 p 2 ,a 2 A 2 ) . The lattice 
momentum can be developed as: 

a 2 p 2 = a 2 /? 2 + c\a 4 p^ -\ (3.6) 

with c\ a constant that depends on the discretization chosen. Then: 

4 / [4] \ 

a 2 p 2 = Y,a 2 pl=a 2 p 2 + c 1 a 4 p [4] + --- = a 2 p 2 \ I +c 1 a 2 ^- T + ■ ■ ■ (3.7) 
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where = YTu.=i P\i- If tne lattice spacing is small, e = rp" /p«\ and we can develop 2 in 
powers of e: 



[4] 



2 A 2 



e(aV,« 2 A 2 ) + 



dQ 
de 



2P [4] , 

a -^- + - 

e=0 P 2 



(3.8) 
(3.9) 



H(4) methods are based on the appearance of a 0'{a 2 ) corrections driven by a term. The 
basic method is to fit between the whole set of orbits sharing the same p 2 the coefficient R and the 
extrapolated value of Q free from H(4) artefacts. In particular we assumed that the coefficient 



dQ (aV(l + cie + ---),« 2 A 2 ^ 



R(a 2 p 2 ,a 2 K 2 ^ 



de 



e=0 



has a smooth dependence on a 2 p 2 over a given momentum window. This can be achieved by 
developing R as R = Rq +R\a 2 p 2 and making a global fit in a momentum window between (p — 
8,p + 8) to extract the extrapolated value of Q for the momentum p and shifting the window for 
every lattice momentum. This procedure of fitting is somehow different from the previous one, 
since the extrapolation does not rely on any particular assumption for the functional form of R. On 
the other, the systematic error coming from the extrapolation can be estimated by modifying the 
width of the fitting window. 



4. Results 



4.1 Calibration of lattice spacings 



The running of aj given by the combination of Green functions in eq. ( |2.4[ ) does depend 
in principle on the momentum and the cut-off. Nevertheless, if we are not far from the continuum 
limit, and discretization errors are treated properly, the coupling will depend only on the momentum 
(except, maybe, finite volume errors at low momenta). 

The procedure to compute the ratio of lattice spacings is then straightforward: it can be ob- 
tained by requiring the estimates of a r for two different simulations (two different j8's) to match 
properly each other. This method has proven to be successful in quenched lattice simulations [1], 
with a deviation with respect to usual Sommer parameter estimates lower than 5%. 

The Nf = 2 results can be seen in figure |4.1| , where the lattice spacing for the lower j8 (j8 = 3.9) 
has been assumed to coincide with the one given in [6] and the other two, for j8 = 4.05 and j8 = 4.20 
are fitted to match the data. 

The deviations are found to be smaller than 5% (see Tab. fP| ), as in the quenched case. This 
deviation might be a signal of discretization errors still present at these j3's. Another source of 
discrepancy could be a possible dependence of results on the quark masses. Further efforts should 
be done in this sence. 
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This paper 


Sommer scale 


deviation (%)) 


a(3.9)/a(4.05) 


1.223(3) 


1.277 


4.2 


a(3.9)/a(4.2) 


1.503(5) 


1.547 


2.9 



Table 1: Best-fit parameters for the ratios of lattice spacings. The error is purely statistics. 




0.3 1 1 1 1 1 1 1 1 1 1 0.22 1 1 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 7 8 9 10 10 12 14 16 18 20 22 24 26 28 30 



Figure 1: QCD coupling defined by from the three lattice data sets employed: red squares stand for j3 = 4.2, 
green ones for = 4.05 and blue for j3 = 3.90. Right (left) plot shows estimates for momenta above (below) 
10 GeV 2 . The physical value (in GeV) of the momentum in x-axis is obtained by applying the ratios of 
lattice sizes in tab.l anda(3.9)" 1 = 2.301GeV. 



4.2 A^and (A 2 ) condensate 

The value of A^ can be obtained by inverting (2^) for the lattice values of aj obtained 



from the lattice for each momentum. When done (figure g) the values of A^ obtained have a 
strong dependence on the momentum, showing the presence of some non-perturbative effects not 



taken into account in (2.5 ). The values of Aj^ are around 320 — 360MeV, much higher than other 
estimations. 

The first non-perturbative correction that does appear un Landau gauge is the (A 2 ) gluon con- 
densate, whose effects on the running coupling are included in ( |2.9| ). The values of A^ and (A 2 ) 
can be simultaneously fixed from lattice data using, for example, the "plateau" method, shown 
in [1]. It consist in varying the value of the condensate to look for a "plateau" in A^ over a given 
momentum window. 

In fig. ||, we also plot A^- derived from confronting the lattice value of aj with the pertur- 
bative+OPE prediction, in terms of the momentum where <Xj is estimated from the lattice. The 
application of the "plateau" method allows us to get as a best estimate: 

A^ = 267±llMeV ; (4.1) 

where again the error takes into account no systematic effect. This result is in good agreement with 
other estimations in litteraure [13]- [15]. The value of the (A 2 ) obtained is 

S 2 r(A\ M =9.6±0.6GeV 2 ( 4 - 2 ) 
which shows a significant increase respect to previous quenched estimates [1]. 
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Figure 2: A-^g derived from fitting the lattice value of (Xt with the perturbative+OPE prediction, in terms of 
the momentum where aj is estimated from the lattice, as described in ref. [1]. 

5. Conclusions and outlooks 

We calculated the running coupling in the Taylor scheme with Nf = 2 flavours of dynamical 
quarks. We found that the matching of the results obtained for different j3's allows to compute the 
ratio of lattice spacings, with a deviation with respect to the string tension always smaller than 5%. 

By comparing the lattice result with the expectation coming from perturbation theory, we 
found the need for a dimension-two gluon condensate associated to a non-perturbative power cor- 
rection. Including this term allows for an agreement between lattice and continuous formulae and 

Nf=2 

then the extraction of the scale . Our result is in agreement with previous determinations. 

The application of this method is straightforward for a higher number of quark flavours and 
might be used in forthcoming Nf = 2 + 1 + 1 lattice simulations. 

As an outlook, we are interested in checking the mass-dependence of our results. In particular 
two effects are to be expected. The first one, at the level of the calibration, could show a dependence 
of the lattice spacing both on j3 and }X. In any case this should not affect our results. The second 
one could be the effect of the mass on the coupling, which seems to be rouled out because of the 
good overlap of the coupling already observed at different /jl's. 
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